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etecting  edges  is  an  important  first  step  in  the  solution  of  many  image  ana¬ 
lysis  tasks.  Edges  are  usedprimarily  to  aid  in  the  segmentation  of  an  image  in¬ 
to  meaningful  regions,  but  are  also  extensively  used  to  computer  relatively 
local  measures  of  textural  variation  (whic,  of  course,  can  subsequently  be  used 
for  segmentation  purposes).  Although  there  has  been  a  considerable  amount  of 
research  concerning  quantitative  models  for  edge  detection,  very  little  work 
has  been  devoted  to  developing  such  models  for  images  described  by  texture 
models . 
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1.  Introduction 


Detecting  edges  is  an  important  first  step  in  the  solution  of  many 
image  analysis  tasks.  Edges  are  used  primarily  to  aid  in  the  segmentation 
of  an  image  into  meaningful  regions,  but  are  also  extensively  used  to 
compute  relatively  local  measures  of  textural  variation  (which,  of  course, 
can  subsequently  be  used  for  segmentation  purposes).  Although  there  has 
been  a  considerable  amount  of  research  concerning  quantitative  models  for 
edge  detection  (e.g.,  Nahi  [1],  Modestino  and  Fries  [21,  Shanmugan  et  al 
[3),  Cooper  and  Elliot  (4]),  very  little  work  has  been  devoted  to 
developing  such  models  for  images  described  by  texture  models.  This  paper 
addresses  the  problem  of  detecting  edges  in  what  are  called  macro- textures, 
i.e.,  cellular  textures  where  the  cells,  or  texture  elements,  are 
relatively  large  (at  least  several  pixels  in  diameter). 

Once  edges  are  detected  in  texture  regions,  they  can  be  used  to 
define  texture  descriptors  in  a  variety  of  ways.  For  example,  one  can 
compute  "edge  per  unit  area"  (Rosenfeld  [5]).  More  generally,  one  can 
compute  first-order  statistics  of  edge  properties  (6,7],  such  as 
orientation,  contrast,  fuzziness,  etc.,  or  higher-order  statistics  which 
can  measure  the  spatial  arrangement  of  edges  in  the  texture.  Such 
statistics  can  be  computed  from  generalized  cooccurrence  matrices  (Davis 
et  al  (8,9))  which  count  the  number  of  times  that  specific  pairs  of  edges 
occur  in  specific  relative  spatial  positions.  Clearly,  the  utility  of 
such  tools  depends  on  the  reliability  with  which  edges  can  be  detected 
in  textures. 

This  paper  is  organized  as  follows:  Section  2  contains  a  description 


of  the  imago  texture  models  which  will  be  considered.  These  models  are 
one-dimensional,  since  the  edge-detection  procedures,  described  in 
Section  3,  are  one-dimensional.  Section  4  contains  derivations  of  the 
expected  value  and  variance  of  the  edge  operator  described  in  Section  3 
and  describes  optimal  edge  detection  procedures  based  on  that  analysis. 
Finally,  Section  5  contains  conclusions  and  a  summary. 
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2.  Texture  Models 


There  are  a  large  number  of  formal  image  texture  models  which  have 
been  proposed  and  studied  during  the  past  few  years.  These  can  be 
broadly  classified  as  pixel-based  and  region-based  models  (Ahuja  [10]). 
All  of  these  models  treat  textures  as  a  two-dimensional  phenomena,  which 
is  appropriate  for  many  applications  (e.g.,  some  medical  applications, 
geographical  applications).  However,  for  other  applications,  regarding 
textures  in  this  way  is  inappropriate;  one  should,  instead,  model  the 
texture  as  a  surface  in  space  with  certain  reflectance  properties.  An 
image  of  such  a  texture  is  then  determined  by  the  spatial  disposition  of 
the  surface  and  the  viewer,  the  frequency  response  of  the  viewer,  and  the 
positions  of  all  light  sources.  Horn  [11]  should  be  consulted  for  an 
introduction  to  this  branch  of  image  science.  Such  models  will  not  be 
considered  in  this  paper. 

,  Pixel-based  models  are  ordinarily  time-series  models  or  random  field 
models.  Time-series  models  have  been  investigated  by  McCormick  and 
Jayaramamurthy  [12],  and  by  Tou  et  al  [13].  Random  field  models  are 
discussed  in  Wong  [14],  and  Pratt  et  al  [15].  For  more  references,  see 
[16-17]. 

This  paper  will  be  concerned  with  region-based  texture  models.  In 
particular,  we  will  consider  one-dimensional  models  which  are  related  to 
two-dimensional  cell  structure  models.  Cell  structure  models  describe 
textures  as  mosaics,  and  can  be  generated  by  the  following  two-step 
process: 

1)  A  planar  region  is  tessellated  into  cells,  ordinarily  convex. 


2)  Each  cell  is  independently  assigned  one  of  m  colors, 
using  a  fixed  set  of  probabilities,  p^,. ■ • »Pm* 

This  process  partitions  the  original  region  into  subregions,  which 
are  the  unions  of  cells  of  constant  color.  If  A  is  the  original  region, 
then  Aj,...,A  are  the  subregions.  Note  that  the  simple  colors  can  be 
replaced  by  more  complex  coloring  processes,  e.g.,  the  grey  levels  in  a 
cell  can  be  chosen  according  to  a  given  distribution,  d,  which  is  itself 
chosen  from  a  set  of  distributions,  D,  according  to  the  given  probability 
vector,  P.  Ahuja  [10]  contains  an  extensive  survey  of  such  models. 

We  will  consider  a  similar  class  of  one-dimensional  models.  A 
texture  model  is  an  ordered  pair  <P,C>  where 

1)  P  is  a  cell  width  model,  which  successively  drops  intervals 
along  a  line,  and 

2)  C  is  a  coloring  model,  consisting  of  coloring  processes, 
cr*,”cm’  and  Probabil ities,  Pi»..**Pm*  As  p  produces  cells,  C  colors 
the  cells. 

If  we  let  w  be  the  random  variable  corresponding  to  cell  width,  then 
the  following  are  examples  of  cell  width  models: 

1)  Constant  cell  width  model 

(  1  w=b 

Pc(w)  = 

c  (  0  w/b 

2)  Uniform  cell  width  model 


3)  Exponential  cell  width  model 
P£(w)  =  A  exp  (  Awl  . 

To  simplify  the  analysis,  we  assume  that  there  are  only  two  coloring 
processes,  c^  and  C2<  Therefore,  there  is  only  one  relevant  probability 
for  chqosing  cell  colors,  which  will  be  denoted  by  p.  Each  coloring 
process  colors  a  cell  by  choosing  the  intensity  of  each  point  in  that  cell 
independently  from  a  normal  distribution  of  intensities.  The  distributions 
are  denoted  N(m.  .v^,  i  =  1 ,2 ,  with  mean  and  variance  v^ 

Notice  that  given  a  one-dimensional  cell  structure  model,  <P,C>,  one 
can  derive  a  one-dimensional  component  structure  model,  where  a  component 
is  a  contiguous  set  of  identically  colored  cells.  For  example,  for  the 
cell  model  <Pc,C>  the  corresponding  component  model  has  components  whose 
lengths  are  distributed  geometrical ly.  Of  course,  in  a  component  model, 
the  various  colors  alternate,  since  by  definition  two  adjacent  components 
must  have  different  colors.  The  component  model  is  required  to  determine 
the  prior  probabilities  of  various  types  of  pixels  (see  Section  4). 


3.  A  One-Dimensional  Edge-Detector 


In  this  section  we  will  describe  a  simple,  one-dimensional  edge 
detection  procedure.  There  are  a  variety  of  reasons  for  considering 
one-dimensional  edge  detectors: 

1)  Computational  efficiency  on  conventional,  sequential  computers. 
Since  we  are  interested  in  detecting  edges  for  the  purpose  of  describing 
texture,  it  is  important  that  the  edge  detection  process  be  made  as 
efficient  as  possible. 

2)  Suitability  for  implementation  on  special-purpose  image 
processing  hardware.  There  has  been  a  significant  amount  of  research 
and  development  of  image  processing  hardware  over  the  past  few  years 
(CLIP  [18],  PICAP  [19]).  One-dimensional  edge  detectors  can  be  easily 
implemented  on,  for  example,  series-parallel  machines,  where  a  row  of 
processors  "scan"  and  process  the  image  one  row  at  a  time. 

3)  Mathematical  tractabil ity.  Direct  analysis  of  two-dimensional 
edge  detectors  is  complicated  by  many  factors,  including  the  mathematical 
complexity  of  most  two-dimensional  cellular  texture  models  and  the 
greater  variety  of  edge-like  features  in  two  dimensions. 

The  class  of  edge  operators  which  we  have  considered  is  based  on 
differences  of  averages  between  adjacent,  symmetric  one-dimensional  image 
neighborhoods.  Specifically,  let  f  be  a  one-dimensional  image  function. 
Then  the  edge  operator  is: 

k 

ek(i)  =  (1/k)  l  (f(i-j)-f(Hj)) 

*  j=l 


=  (1/k)  (LS(I)-RS(I)) 


where 


k 

LS( i )  =  l  f(i-j)  and 
j=l 

k 

RS ( i )  =  l  f(i+j) 
j=l 

By  noticing  that 

ek(i+l)  =  ( 1/k)  [LS(i)-f(i-k)  +  f(i)-RS(i)  +  fl  i+l)-f( i+k+1) ] 
we  see  that  ek  can  be  computed  in  a  constant  number  o1'  operations  per 
picture  point,  independent  of  k,  on  a  conventional  sequential  computer. 

The  operator  ek  is  used  to  detect  edges  by  the  following  three  step 
process: 

1)  Compute  ek(i)  for  all  points  i. 

2)  Discard  all  i  with  | ek ( i ) |  <  t .  This  thresholding  step  is  intended 
to  discriminate  between  points  which  are  edges  of  texture  elements  and 
points  which  are  in  the  interior  of  texture  elements,  but  far  from  edges. 

3)  Discard  all  i  with  |ek(i)|  <  |ek(i+j)|,  |j|  <d.  This  non-maxima 
suppression  step  is  intended  to  discriminate  between  points  which  are 
edges  and  points  which  are  interior  to  texture  elements,  but  are  close 

to  edges. 

Step  3  is  crucial  since  ek  gives  high  response  not  just  at  edges, 
but  also  near  edges,  so  that  thresholding  alone  would  result  in  a  cluster 
of  detections  about  each  true  edge  point.  The  above  procedure  involves 
three  classes  of  texture  pixels: 

1)  edge  pixels,  which  are  pixels  located  directly  at  the  edges 
between  texture  elements; 

2)  near-edge  pixels,  which  are  located  within  distance  d  of  an 


edge  pixel  and  are  discarded  by  the  non-maxima  suppression  step;  and 

3)  interior  pixels,  which  are  located  distances  greater  than  d 
from  the  nearest  edge  pixel,  and  are  ordinarily  eliminated  by  the 
thresholding  step  (but  may  be  eliminated  due  to  proximity  to  above 
threshold,  near-edge  pixels). 

Optimizing  the  above  edge  detection  procedure  involves  choosing  k, 
t  and  d  in  order  to  minimize  the  probability  of  error — i.e.,  minimizing 
the  frequency  of  discarding  edge  points,  and  not  discarding  near-edge 
and  interior  pixels.  This  paper  considers  the  edge-interior  discrimination 
problem  only.  Therefore,  we  will  be  concerned  with  choosing  values  for  k 
and  t  only.  The  complete  edge  detector,  including  the  non-maxima 
suppression  step,  is  discussed  in  Davis  and  Mitiche  [20]. 


4.  Analysis  of  ek 

In  this  section  we  will  derive  the  expected  value  and  variance  of 
ek  at  edges  and  at  interior  points.  We  will  regard  an  interior  point 
as  a  point  whose  distance  from  the  nearest  edge  is  greater  than  k. 
Expressions  for  the  prior  probabilities  of  edge  and  interior  points  are 
developed.  Finally,  by  assuming  that  e^  is  normally  distributed  at 
edges  and  interiors,  a  minimum  error  thresholding  procedure  for 
distinguishing  between  edge  points  and  interior  points  is  developed. 

4.1  The  expected  value  of  ek,  E [ ek ] 

The  definition  of  ek  was  originally  given  for  a  discrete  function  f. 
If  f  is  continuous,  then  we  can  redefine  ek  as 

(0  i  k 

ek(i)  =  1/k  f(i+j)dj  -  I  f(i+.j)dj 
K  -k  0 

Then,  the  expected  value  of  ek  is 

r  r°  rk 

E[e.(i)J  =  1/k  E[f(i+j)Jdj-  E[f(i+j)]dj 
K  l  -k  0 

If  i  is  an  interior  point,  then  all  points  f(i+j),  -k<j£k  are  colored 
by  the  same  process.  Therefore,  the  expected  values  are  all  the  same, 
and  thus  E f ek ( i ) J i  is  an  interior  point]  =  0. 

Now  suppose  that  i  is  an  edge  point,  and  assume,  without  loss  of 
generality,  that  the  cell  to  the  left  of  i,  C^,  is  colored  by  process  c^, 
and  that  the  cell  to  the  right  of  i,  Cr>  is  colored  by  by  process  C£»  and 
that  mj  xr^.  Let  w ^  be  the  width  of  and  wr  be  the  width  of  Cr  (see 
Figure  1).  We  will  also  make  the  simplifying  assumption  that  the  points 


to  the  left  of  C(,  (or  the  riqht  of  f  )  are  individually  colored  with 
processes  c^  and  c 2  with  probabilities  p  and  (1-p).  For  k  much  greater 
than  or  wr,  this  assumption  is  not  unreasonable.  As  or  wr 
approaches  k,  it  is  more  likely  that  only  one  cell  will  be  found  to  the 
left  of  or  the  right  of  Cr.  However,  large  cells  are  ordinarily  less 

likely  than  small  cells.  Letting  a  =  prn^  +  ( l-p)^ ,  we  can  then  write 

E  [ek(  i )  ]  =  1/k  ([  (m1w£+a(k-w^)P(w^)dw^  +  J  m1kP(w^)dw£) 

f  (m2wr+a(k-wr)P(wr)dwr +  J  m2kP(wr)dwr) 

0  k 

Since  w^  and  w^  are  drawn  from  the  same  distribution,  terms  can  be  grouped 
to  obtain 

r  fk  r°°  -1 

F[ek(i)]  =  1/k  J  (m1-m2)wP(w)dw + j  (n^-n^) kP(w)dw 

r  r00  r00  i 

=  1/k  I  (m1-m2)w0- J  (m1-m2)wP(w)dw +  J  (m1-m2)kP(w)dw 
k  k 

(mrm?)  r  foo 

=  -  w  (w-k)P(w)dw 

k  L  u  J.  j 


where 


wn  =  wP(w)dw. 

u  JQ 


4.2  The  variance  of  e, 

k 

From  the  definition  of  e,  ,  we  have 

k 


T  fk  1 

Var  [e. ( i ) J  =  Var  1/k  f ( i-j)-f( i+j)dj 

*  L  n  J 


Suppose  i  is  an  interior  point  and  is  in  a  cell  colored  by  c^.  Then 

Var[ek(i) ]  =  (1/k2)  k  2vj 
=  2vj/k 

If  i  is  an  interior  point  and  is  in  a  cell  colored  by  c2,  then 
Var[ek( i ) ]  =  2v2/k. 

Next  we  will  consider  the  case  of  i,  an  edge  point.  Let 
LS(  i )  = 


k 

f(i+j)dj 


RS(i)  =  J 


Then 

ek(i)  =  1/k  (LS(i)-RS(i)) 

We  will  derive  expressions  for  Var [ LS ( i ) ]  and  Var [ RS ( i ) ) .  There  are  two 
cases  to  consider: 

1)  w^  or  wf  >  k 

2)  w^  or  wr  <  k 

If  w^k,  then  the  variance  of  LS ( i )  is  kv^.  If  wr>.k,  then  the  variance 
of  RS ( i )  is  kv2- 

Next,  suppose  w^<k.  We  will  assume  that  of  the  k-w^  pixels  not 
in  C^,  i j  pixels  are  colored  by  process  Cj  and  £2  =  (k-w^)-£^  pixels  are 
colored  by  process  c2,  where  is  a  random  variable  described  by  a 
binomial  distribution  with  parameter  p  and  f2  is  a  random  variable 
described  by  a  binomial  distribution  with  parameter  (1-p). 

In  general,  if  y  is  the  sum  of  a  random  number,  n,  of  independent 


experimental  values  of  a  random  variable  x,  then 
Var[y]  =  E[nJ  Var[x]  +  (E[x] )2  Var[n] 

Thus,  if  w^  is  a  fixed  value  less  than  k, 

Var(LS(i) |w£ <  k]  = 

w£vl  +  El ^1 1 V1  +ml  Varl'ti1 
+  E(f2]v2  +  m2  Var[£2 ] 

with 

E [£x 1  =  (k-w£)p 

E[£2)  =  (k-w£)(l-p) 

Var[£1l  =  Var[f2]  =  (k-w£)p(l-p) 

A  similar  expression  can  be  obtained  for  RS(i). 

Combining  the  cases  w^k  and  w^<k,  we  can  write 

Var[LS( i ) | i  an  edge]  = 

Tz[[  kvlP(w£)dwC 
k  L  k 

fk 

+  J  (w£v1  +  (k-wc)(pv1  +  {l-p)v2) 

+  (k-w£)  p(  1-p) (m2  +  m| ) )  P(wf)dwf 

Again,  a  similar  expression  is  obtained  for  RS ( i ) .  Since 
Var[ek( i ) ]  =  1/k2  [ Var[lS( i ) ]  + Var[RS( i ) ] J 


Since  w^  and  wr  are  drawn  from  the  same  distribution,  we  group  terms  and 


m  fmmmm mm. 


pg,  and  of  an  interior  point,  p..  An  edge  point  is  defined  to  be  a  point 
on  the  image  line  which  is  no  more  than  a  fixed  distance  d  away  from  a 
true  edge.  We  choose  d>0  so  that  pi  to.  An  interior  point,  as  before, 
is  a  point  that  is  at  least  distance  k  away  from  a  true  edge. 

In  the  following  we  will  consider,  without  loss  of  generality,  only 
points  within  Cj-components .  Given  that  a  point  is  picked  at  random, 
let  h  be  the  probability  density  function  for  the  distance  y^  to  the  next 
edge  (the  forward  edge  or  edge  on  the  right).  From  general  results  on 
random  incidence  into  a  renewal  process  we  have  [21] 

h(y  j )  =  ( 1  -  p  r  o  b[r<^y1])/E[r] 

The  function  h  is  also  the  probability  density  function  for  the  distance 
y2  to  the  backward  edge  (the  nearest  edge  on  the  left  of  the  point). 

Then  we  have 


hj(yry2)  =  ^(y1)hc(y2/y1) 

where  h  is  the  conditional  probability  function  for  y„  and  h.  is  the 
c  c  j 

joint  probability  function  for  y^  and  y^.  If  y^  and  y^  are  independent, 
then 

Vyry2^  =  My^Myg) 

From  the  above  we  now  can  derive  expressions  for  the  priors 


For  the  discrete  case  we  have 


P-j  =  I  I  h(y1,y2) 
1  yj-k  y2=k  1 


Pe 


-  1-  l 

yx=d+l 


l  h(y,,y2) 

y2=d+l 


The  set  of  interior  points  and  the  set  of  edge  points  are  not  complementary 
in  the  sense  Jhat  they  do  not  form  a  partition  of  the  set  of  all  points  of 
the  image  line.  Thus  it  is  necessary  that  we  normalize  pe  and  p... 

pe  *  pe/(petpf) 

p,  *  lyip.+p,) 


Examples 

1.  Geometric  distribution  of  component  widths.  Assume  for  simplicity 
that  b  =  1.  Then 

yl 

h(y,)  =  (1-P)(1-  l  (l-p)py_1) 

1  y=i 

=  U-p)  [  ( i-p)py_1 

y>yx 

=  ( l  -  p  )  pY  1 

co  y, 

Prob[y,  >  k]  =  l  (l-p)p 

yfk 

oo  y  .-1 

=  p  l  (l-p)p  1 

yfk 


Thus 


Elr]  =  E[n]  E[w] 


Var[r]  =  E[n]  Var[w]  +  (E[wJ )2  Var[n J 

If  w  Is  continuous  and  if  f  and  g  are  the  probability  density  functions 
for  w  and  r,  then 

gT(s)  =  P /  (fT(s)) 

where  f^  and  are  the  exponential  transforms  (s-transforms)  of  f  and  g 
respectively,  and  p^  is  the  z-transform  (discrete  transform)  of  p^. 

If  w  is  discrete,  then  the  above  expression  holds  with  f  and  g 
being  the  probability  mass  functions  for  w  and  r,  and  fT  and  gT  being 
their  respective  z-transforms.  By  taking  the  inverse  transforms  one  can 
obtain  the  distribution  of  component  widths. 

Examples 

1.  Constant  distribution  of  cell  width.  For  the  constant  cell  width 
model,  g  is  described  by  the  geometric  probability  mass  function 

/  £  - 1 

l(l-p)p°  for  r  =  b,2b,3b,... 

9(r)  =  j 

(  0  otherwise 

E [ r ]  =  b/( 1-p) 


Var[r]  =  b2  p/( 1-p) 2 


2.  Exponential  distribution  of  cell  width.  For  the  exponential  model 
with  parameter  A,  g  is  described  by  the  probability  density  function 


9(r)  =  (1-P)  A  Exp(-(l-p)  Ar) 

Thus  component  widths  are  still  exponentially  distributed  and 
E[r)  =  l/( ( 1-p) A ) 


3.  Uniform  distribution  of  cell  width. 


E[rj  =  b/(2( 1-p) ) 


Var[r  1  = 


b/(4( 1-p) )  _ 

( 1/3  +  p/( 1-p) ) 


and 

9(r)  ■  e(p/b)rg'(r) 


where 


g'(r) 


(1-P) 


y  (-■)Ve"IP  (r-bn)" 
n--0  bntl  iw 


bn 


(r) 


i 

n=0 


,  , x n  n  -np 
(-D  P  e 
i.n+1 


[r-b(n-t-l)  ]n 

r(n+l) 


Sb(n*l)(r) 


and 


sk(t)  = 


1 0 
I* 


0  <  t  <  k 
t  >  k 


is  the  unit  step  function. 


In  order  to  perform  minimum  error  thresholding  on  a  given  image 
texture,  it  is  necessary  to  know  the  prior  probabilities  of  an  edge  point. 


finally  obtain,  for  edge  points, 

Var(ek(i) | i  is  an  edge]  = 

1/k^  j^k(vj+v2)  J  P(w)dw 

r  k 

+  (v,+v?)  wP(w)dw 

1  c  0 

fk 

+  2(pv,+(l-p)v,)  (k-w)P(w)dw 
1  c  0 

+  2p(l-p)(m^ +m2)  j  (k-w)P(w)dwj 

For  tl.e  exponential  model  with  p=  .5,  for  example,  we  find  that 

(m  -m  )(l-e'Xk) 

Elek(i)li  is  an  edge)  =  - - 

and 

i  ^vl+v2) 

Var[e.(i)|i  is  an  edge]  =  -  + 

K  k 

Other  simple  expressions  can  be  obtained  for  the  constant  and  uniform 
models. 

In  order  to  derive  a  minimum  error  threshold  for  discriminating 
between  edge  and  interior  points  using  e^,  it  is  necessary  to: 

1)  determine  the  prior  probabilities  of  edge  and  interior  points, 
and 

2)  specify  a  form  for  the  distribution  of  e^  at  edges  and  at 


2  ,  2 
m^  +  m2 

2k  _ 


(> 


1-e 


-Ak 


Ak 


interior  points. 


To  compute  the  prior  probabilities  of  edge  and  interior  point,  we 
must  derive  a  component  model  from  the  cell  width  model.  A  component  • 
is  a  set  of  connected,  identically  colored  cells.  Let  a  Cj-component 
(Cj-component)  be  a  component  whose  cells  are  colored  by  process  Cj  (C2). 
Then,  the  length  of  a  component  is  defined  to  be  the  number  of  cells 
that  compose  it,  and  the  width  of  a  component  is  its  actual  measure. 

In  the  following  we  will  show  how  a  component  model  can  be  derived 
from  a  cell  model.  As  mentioned  before,  the  cell  coloring  process  is 
a  Bernoulli  process  which  selects  coloring  process  c^  with  probability  p 
and  coloring  process  c^  with  probability  ( 1 - p ) .  Therefore,  the 
c ^-component  and  C2-component  lengths  are  random  variables  described  by 
a  geometric  probability  mass  function.  In  the  following  analysis  we  will 
consider  only  c ^-components  and  simply  refer  to  them  as  components. 

The  same  analysis  will  hold  for  C2~components. 

If  n  is  the  random  variable  that  represents  component  length  and 
p^(n)  is  the  probability  mass  function  for  it,  then 

P^(n)  =  (I-p)  pn_1 

Now,  if  w  and  r  are  random  variables  that  describe  the  cell  width 
and  the  component  width  respectively,  then  r  is  the  sum  of  a  random 
number  n  of  independent  identically  distributed  experimental  values  of  x: 
n 

r  "  l  *. 

i  =  l  1 

Expressions  for  the  expected  value  and  variance  of  r  can  be  obtained  in 
terms  of  the  expected  value  and  variance  of  n  and  w: 


Since  the  geometric  distribution  is  memoryless,  y^  and  y2  are  independent 
random  variables,  and 

Pj  =  Probly^k,  y2>_k]  =  pk  pk  =  p2k 
p,  - 1  -  p2(dtl) 

Normal izing. 


p,  =  P2k/(p2k*i-P2(d+1)> 


pe  •  i-P, 


2.  Exponential  distribution  of  component  widths.  Let  \  be  the  parameter 
of  the  distribution.  Then 

1  -  (1  -  Exp(-Ay. )) 

h(yL)  =  - — - i—  =  \  Exp(-*yi) 

Prob[y^k]  =  Exp(-Xk) 

The  exponential  distribution  is  also  memoryless  so  that  y^  and  y2  are 
independent.  We  then  have 

Pi  =  Exp(-2Ak) 

P„  =  1  -  Exp(-2d\) 


Normalizing, 


P1  =  Exp(-2Xk)/(Exp(-2Xk)  +  1  -  Exp(-2dX) ) 


=  1  -  P- 
^e  ri 


In  order  to  use  pg,  p.  and  the  expected  values  obtained  above  for 
minimum  error  edge  detection,  we  will  make  the  assumption  that  ek  is 
normally  distributed  at  edge  points  as  well  as  interior  points.  More 
precisely,  we  assume  that  ek  is  N(0,2v/k)  at  interior  points  and 
N{E[ek),Varfek])  at  edges.  This  assumption  is  certainly  valid  for 
interior  points  since,  in  this  case,  LS  and  RS  are  each  the  sum  of  k 
independent  experimental  values  sampled  from  the  same  normal  distribution. 
At  edges,' each  of  LS  and  RS  is  the  sum  of  k^  independent  experimental 
values  drawn  from  the  normal  distribution  that  represents  one  of  the 
coloring  processes  and  k2  independent  experimental  values  from  the  normal 
distribution  that  describes  the  other  coloring  process.  The  sum  of  k^ 
and  k^  is  k,  but  kj  and  will  in  general  vary  from  point  to  point  on 
the  image  line. 

The  assumption  of  normality  of  ek  at  edges  will  hold  well  if  the 
variances  of  k ^  and  k^  are  low.  This  means  that  in  the  neighborhood  of 
any  edge  point  on  the  image  line  the  number  of  pixels  colored  by  either 
process  remains  almost  constant.  For  example,  this  is  trivially  true  for 
the  constant  distribution  of  cell  widths,  with  k£b,  b  being  the  width  of 
a  cell.  In  this  extreme  case,  kj  =  k  and  k^  = 0  or  vice  versa.  However, 
if  kj  and  k£  have  a  high  variance,  then  the  normality  assumption  will  not 
hold  very  wel 1 . 

Thus,  an  important  factor  for  the  validity  of  the  assumption  is  the 
variance  of  cell  widths.  Ideally,  this  variance  should  be  small;  however 
another  property  that  would  tend  to  make  the  assumption  hold  well  is  that 
the  image  model  be  more  likely  to  contain  cells  whose  widths  are  close  to 
k.  This  will  keep  the  probability  that  the  neighborhoods  and  Cp  extend 


over  more  than  one  cell  quite  low. 

Given  the  normality  assumption,  the  following  two  step  process  can 
be  used  to  compute  an  optimal  k  and  t  for  a  minimum  error  edge  detector. 

1)  For  a  range  of  k,  find  the  minimum  error  threshold  for 
discriminating  between  edges  and  interior  points.  Since  both  ek  at 
edges  and  ek  at  interior  points  are  modeled  by  normal  distributions  with 
known  parameters  and  priors,  this  is  straightforward.  Let  er(k)  be  the 
probability  of  error  for  the  minimum  error  threshold  for  e^  and  let  t(k) 
be  the  threshold. 

2)  Choose  k  such  that  er(k)  <_er(k'),  for  all  k'  considered.  Then 
(k,t(k))  define  the  minimum  error  edge  detector. 

Figure  2  shows  plots  of  er(k)  as  a  function  of  k  for  the  three  cell 
width  models  presented  in  Section  2.  Notice  that  the  value  of  k  which 
minimizes  total  errors  is  the  mean  cell  width  for  all  three  models.  The 
reason  that  the  curves  tend  to  level  off  at  high  k,  rather  than  rise  to 
a  high  error,  is  that  as  k  becomes  very  large,  the  prior  probability  of 
interior  points,  pi ,  approaches  zero.  For  very  high  values  of  k,  the 
near  zero  value  of  p,.  causes  the  programs  which  compute  minimum  error 
thresholds  to  become  unstable.  Therefore,  we  arbitrarily  stopped 
computing  er(k)  for  k > 20.  Figure  3a-e  shows  an  example  of  the  effect 
of  k  on  the  performance  of  e^.  Figure  3a  contains  a  checkerboard  texture 
with  b=16,  p=.5,  m1=30,  m^  =  20 ,  and  v1=v2=10.  Figure  3b  shows  the 
true  edges,  while  Figures  3c-e  show  the  results  of  applying  e^,  k  =  8,  16 
and  24,  thresholding  at  the  minimum  error  threshold  for  the  appropriate  k, 
and  then  performing  non-maxima  suppression  across  8  pixels.  Note  that  the 
results  for  the  optimal  value  of  k  (16)  are  significantly  better  than 
choosing  k  too  small  (8)  or  too  large  (24). 


5.  Discussion 


This  paper  has  discussed  the  problem  of  detecting  edges  in  cellular 
textures.  A  general  edge  detection  procedure  was  proposed.  This  procedure 
involved  applying  an  edge  sensitive  operator  to  the  texture,  and  then 
thresholding  the  results  of  the  edge  operator  and  finally  computing  "peaks" 
from  the  above  threshold  points.  This  paper  concentrated  on  the  thresholding 
process  and  developed  a  minimum  error  thresholding  procedure  based  on  an 
analysis  of  the  edge  operator  e^.  The  thresholding  procedure  assumed 
that  e^  was  normally  distributed  at  edges  and  at  interior  points. 

The  peak  selection  step  was  not  considered  in  this  paper.  It  is 
discussed  in  [20],  which  also  includes  examples  of  choosing  optimal  edge 
detectors  for  real  textures,  and  a  comparative  classification  study  using 
optimal  and  suboptima  1  edge  detectors. 
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PROBAE 


Figure  2.  Plots  of  er(k)  as  a  function  of  k  for  the  three  cell  width 
models.  All  parameters  are  the  same  for  all  three  models: 
=  10,  =  5 ,  v^=v,,=2,  P  =  . 5  and  Wg  *  10. 


Figure  3.  Application  of  ek  to  a  checkerboard  texture. 


